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Making use of the nonrelativistic potential model for the description of mesons, and working in 
the Shifman-Voloshin limit, we compare the integrated rate T{B — > X c lu) calculated as a sum of 
the individual decay rates to the quantum-mechanical analog of the OPE. In the case of a potential 
regular at the origin, we find a well-defined duality violation, which is however exponentially small. 
It corresponds to the charm resonances kinematically forbidden in the decay process, but apparently 
picked up by the OPE. For singular potentials, we do not obtain a full OPE series, but only a limited 
Taylor expansion, since the cofficients become infinite beyond some order. In this case, we do not 
find an indication of duality violation: the difference is smaller than the last term of the limited 
expansion. This emphasizes that the case of singular potentials, which may be relevant for QCD, 
deserves further study. 

PACS number(s): 13.20 He, 12.39 Jh, 12.39.Pn 



The theoretical framework based on the Operator 
Product Expansion (OPE) determines in QCD the heavy 
meson inclusive decay rate as series in inverse powers of 
the heavy quark mass, with the coefficients proportional 
to the meson matrix elements of the local operators of 
increasing dimensions G|,0|. The calculation is based on 
representing the decay rate as the contour integral in the 
complex go-plane. The OPE makes the contour integrals 
easily calculable term by term and provides the decay 
rate as a 1/toq series. 

There are however potentially dangerous points in this 
calculation: 

(i) the OPE series is at best assymptotically convergent 
even for large absolute values of the complex qo, 

(ii) the integration contour for the decay rate contains a 
segment near the physical region, where the OPE cannot 
be justified Q]. 

This might lead to the violation of duality for the decay 
rate, i.e. to the difference between the OPE-calculated 
decay rate and the result of summing the individual de- 
cay rates of the opened channels. This issue was also 
discussed by N. Isgur || . 

In this letter we discuss the semileptonic decay rate 
in the small velocity limit and use the nonrelativistic 
potential model for the the description of mesons. We 
perform a short-time expansion in operators of increas- 
ing dimensions which we call OPE and which has indeed 
some common features (but also important differences) 
with the OPE expansion in the field theory. We consider 



the two cases: regular confining potentials [] and singular 
potentials. 

In the SV limit both the amplitude and the decay rate 
can be formally obtained as a double expansion in 1 /m c 
and I /5m. We consider lowest orders in l/m c , up to 
l/m\, and all orders in 1/Sm. Note that this involves 
terms of much higher order than usually done when one 
expands in 1/mg with mb/m c , or as well Sm/mQ, fixed. 
Our double expansion allows on the contrary to go much 
further in 1/Sm, and this might allow to display subtle 
duality violations. 

For the regular potential we obtain the full 1/Sm ex- 
pansion, which is only asymptotic to the physical width 
expanded to the same order in l/m c . The difference 
is of order Sm/m^.exp(—Sm/A), which means exponen- 
tially small duality violation. 

For the singular potential we do not obtain the full 
1/Sm expansion: following the same procedure as for 
the regular potential leads to infinite coefficients beyond 
some order in 1/Sm. In this case, we find that the trun- 
cated expansion satisfies duality up to this order. 

We consider the inclusive semileptonic decay B — > 
X c lv in the Shifman-Voloshin (SV) limit A <C 5m = 
mb — m c <C m Cl mf, and treat mesons as nonrelativistic 
bound states of spinless quarks in a confining potential 
(a detailed calculation is given in Q|). This model max- 
imally simplifies both constructing the OPE series and 
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1 A regular potential is a potential which is an analytic func- 
tion of r at r = 0. For example, the potential V(r) ~ \ f\ falls 
out of this class. 

2 As we shall see this expansion contains only a finite number 
of non zero terms. 
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calculating the sum of the exclusive channels. For the 
sake of argument we consider the case of leptons coupled 
to hadrons through the scalar current. In this case the 
leptonic tensor is reduced to a scalar function L(q 2 ). The 
amplitude T depends on the two variables, and we choose 
them as go and q 2 in the S-rest frame: 

T(qo,q 2 ) = \J dxe^(B\T(J(x),J+(0))\B) 

= v \{B\J\X{-q))\ 2 m 
Z-* M B - E x {-q) - qo 1 > 

The sum in (Q) runs over all hadron states with the ap- 
propriate quantum numbers. The states are normalized 
as follows {p\jf} = (2tt) 3 S(p — p 1 ), and E x (-q) is the 
energy of the state X with the total 3-momentum — q. 

At fixed g 2 , T(qo, g 2 ) has a cut in the complex go-plane 
along the real axis for q < Mb — Mb — q 2 /2(m c + rrid), 
see Fig 1. 



A. THE MODEL 

Let us proceed along the lines of ref. Q . We treat the 
leptonic part relativistically, but for the description of 
mesons as bound states of spinless quarks use the nonrel- 
ativistic potential model with a confining potential. We 
consider the decay in the -B-rest frame. The Hamiltonian 
of the bq system at rest has the form 
Hbd = mb+rrid + hbd, hbd = k 2 /2mb + k 2 /2md + Vbd(r), 
such that (h bd ~ £b)\B) = 0, (H bd - M B )\B) = 0, and 
M B = m b + m d + e B - 

The Hamiltonian of the eg system produced in the 
scmilcptonic b — > civ decay reads 

Hcd(q) = m c + m d + (k + q) 2 /2m c + k 2 /2m d + V cd (r). 
The eigenstates of this hamiltonian are \D n (q)) such that 
{H cd (q) - E Dn (q))\Dn(q)) = 0, where E Dn {q) = M Dn + 
q 2 /2(m c + md) and Md„ = m c + m,d + eu„ ■ 
The Qq potential can be expanded as follows: 
V Qq = V Q + Vi/2m Q + V 2 /2m 2 + ... 




a. 

FIG. 1. Singularities of the amplitude T(qo,q 2 ) in 
the complex go-plane. Circles are poles, corresponding to 
low-lying charm states, and the cross marks the location of 
the pole in the free b —> c quark transition. 



A part of this cut for |g| < g < Mb — Md — q 2 /2(m c + 
md) corresponds to the decay process. The decay rate 
can be represented as the contour integral in the complex 
g -plane over the contour C(g 2 ) (Fig 1) 

T(B - X c l V ) = [ dq- 2 \q] [ ^L(q 2 )T{q ,q 2 ). (2) 

J JC(q 2 ) 2m 

The contour C((f) selects at any given q 2 only states 
kinematically allowed in the decay B — » X c lv. It is 
tightly attached to the points P± with the coordinates 
(|g|,±i0), otherwise it can be freely deformed in the re- 
gion where the function To(qo,q^) is analytic. 
The amplitude can be expanded in a series 

TXgcg' 2 ) =^c 4 (g ,g' 2 )( J B|a 4 |B>, (3) 

i 

where Oi are operators of increasing dimensions and 
Ci(go,g 2 ) are the c-number coefficients. Introducing the 
expansion (0) into (||) gives the integrated rate as an OPE 
series^] 



3 The OPE series in the potential model has an important 



B. SUM RULES 

The relationship between the sum over the individual 
channels and the meson matrix elements of the operators 
is established by the sum rules. Let us introduce S n (q) 
through the relation (5 n (q) = £D n —£B — <f'md/2m c (m c + 
md)) 

M B - go - E Dn (q) = Sm-q - q 2 /2m c - 5 n (q). (4) 

The S n (q) is the eigenvalue of the operator 8H(q) 

M B - H cd {q) - 5m - <f/2m c - SH(q) (5) 

with \D n (q)) the corresponding eigenstates. The sum 
rules are obtained by inserting the full system of the 
eigenstates \D n (q)) into (B\ (5H(q)) 1 \B): 

OO 

(B\ (SH(q)y \B) = J2 \FrM\ 2 ■ (6) 

n=0 



distinctions from the Wilsonian scheme in the field theory 
where contributions of distances below the scale are re- 
ferred to the Wilson coefficients, while contributions of dis- 
tances above this scale are referred to the matrix elements. 
As a result, both the Wilson coefficients and the matrix el- 
ements acquire the ^-dependence. In the potential model, 
the T-product of the two currents is also expanded in a se- 
ries of operators of increasing dimensions, but the resulting 
c-number coefficents and the matrix elements in Eq. (Q) are 
scale independent. 
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where F n (q) — (B\D n (q)) is the B — > £>„ transition form 
factor. This relation represents the sum over all cd res- 
onances in terms of the J3-meson matrix element of the 
operators (5H(q)) 1 . For the potential regular at the ori- 
gin r = the sum over n is convergent for any i, whereas 
for the singular potential both sides of Eq. (^) are conver- 
gent for small i and diverge for large i. At the moment 
we proceed formally and discuss this problem in more 
detail in section E. 

C. DUALITY RELATION FOR THE AMPLITUDE 



we can obtain either the width as a sum over the exclu- 
sive final states, or as the OPE series. The expression 
(^) is an accurate approximation to (Jfy only when qo 
is far from the singularities of T(qo,q). The contour C 
can be deformed away from the singularities except near 
its fixed end points. When integrating over qo this is a 
possible source of discrepancy, i.e. of duality violation. 
Consequently, we are now going to estimate the integral 
of expression (Q), i.e. the sum over the exclusive chan- 
nels, and the integral of expression (^), i.e. the OPE 
prediction, and compare both results. 



Making use of the sum rules (jg), we represent the am- 
plitude as a sum of the operators: 



\F n (q)\" 



M B - qo - E n (q) 



(7) 
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(B\(5H(py\)B 
5m ~ £z - lo (5m - £- - qoY 
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(9) 



This expression is the duality relation for the amplitude: 
the sum (^) runs over the infinite number of the charm 
resonances, and the sum (^) runs over the infinite number 
of the operators of the increasing dimensions (the OPE 
series). In fact, the location of singularities in the com- 
plex go-plane in the series (0) and (^|) is quite different: 
in (^) it is an infinite set of single poles at the different lo- 
cations corresponding to different charm resonances, and 
in (JsJ> it is an infinite set of poles of the increasing order 
at the same point. 

However this set of equations is only a formal one ; in 
fact, (0) is a summable series leading to a finite result 
in all cases ; on the other hand, the situation of eq (|^) 
is more subtle. In the singular case, the coefficients are 
infinite beyond some order, and one must accordingly 
truncate the series. In the regular case, the eq (||) is only 
an asymptotic series : notice that the geometric sum over 
i in eq (|^) has a domain of convergence which is repelled 
to infinity with n. 

Let us illustrate it with a simple example: Assume that 
F 2 ~ e~ n and E n ~ n. Then the analog of the above 
equations takes the form 



D. THE OPE CALCULATION OF THE DECAY 
RATE 

Let us first proceed with the amplitude in the form (|^) 
and obtain the OPE expression for the decay rate. We 
consider the leptonic tensor of the general form L(q 2 ) = 
(q 2 ) N . For technical reasons, it is convenient to isolate 
hbd in the expression for 5H(q) as follows 



5H(q) =h bd -€ B + — +[ 

m c \m c m b I 2 



(11) 



Substituting (11) in (|9|) and performing the necessary 
integrations gives a series in 1/ m c |i] 



r OPE (B^ x c iv) 
r(6 -> civ) 



(B\k 2 \B) 



2m 2 . 



-(27V + 3) 



(B\Vi\B) 
2ml 



2AT+3 (~ii+2 

(-1) 4 2Ar + 5 



{B\O t \B) 



i=l 



2N + 5 m 2 r 5m % 



+o 



(12) 



with C l n = ^(^Ij); and O; — k(hhd~ efl)*fc- An important 
feature of the OPE series ( |l2| ) is that the leading-order 
term reproduces the free-quark decay rate, and the first 
correction emerges only in the 1/m 2 order (cf. [0J|]). 



E. SUMMATION OF THE EXCLUSIVE 
CHANNELS 



E 



z — n '■ — ' z 

n=0 n=0 i=0 



Ee ^ f n Y 1 % - 
y ' 4 V y ) y ' 1 yi 



(10) 



The last step is obtained by changing the order of summa- 
tion and using the relation ^ 



n=0 ' 



i* ~ i\ The series 

( |l0| ) in i is only asymptotic and not even Borel summable. 

From the amplitude T under the form eq. ([?]) or eq. 
(^), respectively, by integration over the same contour C, 



Now let us sum the rates of the exclusive channels. 
The B — » D n transition form factors have the form Q 

F 2 {q) = 1 - p\t IK + 0{(?/m 4 c ) + 0{5m 2 l3 2 /m 4 c ), 
Fl% = Plf/ml + Otf/mt) + 0(5m 2 f3 2 /mi). 

Since \q\ < 5m in the decay region, these expressions 
allow calculating the decay rate to the accuracy 5m 2 /m 2 . 
Explicitly, we obtain 0: 
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T(B -> D lv) 



= 1 - 



3p 2 6m 2 



T(b -> civ) 

- (2N + 3) 

T(B -> D w /i/) 3p* Sm 2 
T(b -> di/) 



5m 

2./V + 5 2 l + md/m c to;? 

< B|fc 2 + VMS > 



2m 2 



6m 



2N + 5 w 2 

27V+5 



3 {p 2 n A n ) — 



m 



2 \ .-> y^P ^A^) 



i=2 



2N + 5 6m> 



where A„ = en re - cd ■ 

The main contribution is given by the B —> Dq tran- 
sition. Excited states contribute only starting from the 
(6m) 2 / m 2 order in the SV limit. Notice that each of the 
exclusive rates contains terms of the order Sm 2 /m 2 and 
A6m/m 2 which are absent in the OPE series. 

Summing over all opened exclusive channels gives 



T(B -» XJv) 



T(b -> civ) 



+3 




. Sm ( 

m 2 \1 + md/m c 

B\k 2 + V!\B) 



-(2N + 3) x 



2m 2 



+ 



2N+5 

E 

i=2 



1 \ir<i ^ ( X/ Pn^n 

~ L ) °2AT+5 \ n=l 



2A^ + 5 



m 2 8m l 



(14) 



The sum over the charm resonances is truncated at n max , 
which is the total number of the resonance levels opened 
at q 2 = 0. For the confining potential and in the SV limit 
n max is found from the relation A„ max ~ 5m. 



F. CHECK OF DUALITY FOR REGULAR 
POTENTIALS 

The transition radii in the expression ( |l4| ) are not inde- 
pendent and related to each other through the sum rules. 
These sum rules can be obtained from (^|). Expanding 
both sides of (|6|) in powers of 1 /ttiq and taking the linear 
q 2 term gives the set of the sum rules 0: For i = one 
finds the Bjorkcn sum rule ||, for i = 1 - the Voloshin 
sum rule Q, for i > 2 - higher moment sum rules: 



Using these relations to rewrite the OPE result ( p_2[ ) as 
the sum over hadronic resonances, the difference between 
the OPE and the exclusive sum (the duality-violating 
contribution) explicitly reads 

_ T OPE (B -> X c lv) - T(B -f X c lv) 



T(b -> civ) 



Sm 2 2 ^ 5 (-1)^+5 ^ 2 



3^r- x 

m 



-O 



E 

i=0 

A 2 Sm 



n>n maI 



nr. 



6m 2 3 
^2~/V + 5 



E pi 1 



n>n maa; 



A~<5m 



2JV+5 



(16) 



Quite remarkably, (5r happens to be equal to the sum 
of the extrapolated widths for charm states beyond the 
kinematical limit. A similar expression is found in QCD2 
H . Clearly, the duality- violating effect is connected with 
the charm states forbidden kinematically in the decay 
process. Notice that T OPE (B -> XJv)-T(B -> XJv) < 
0, because A„ > Sm for n > n max , and 2 AT + 5 is odd. 

To estimate the size of the duality- violation effects, the 
behavior of the transition radii and the relation between 
A„ and n max , which will be given by the behaviour of the 
excitation energies also at large n, are needed. For quite 
a general form of the confining potential we can write 
the following relations for A„ at large n A n > ACn a for 



n > n r 



and A r . 



AC(n r , 



Sm, with C and 



a some positive numbers. In particular, this estimate is 
valid for the confining potentials with a power behavior 
at large r. This estimate for A„ is only depending on the 
behaviour of the potential at large distances. 

The behavior of the radii p 2 n at large n are then con- 
nected with the finiteness of the r.h.s. of the sum rules 
(|l5|): for a potential regular at r = 0, the matrix elements 
in the r.h.s. of the sum rules are finite for any i, which 
means that the radii p 2 n are decreasing with n faster than 
any power. Essentially this means that p 2 n ~ exp(— n), 
and therefore the duality- violating effect in the decay rate 
in (|l^) is of order Sy ~ Sm 2 jm 2 c exp(— 6m/ A). One of 
such examples, the harmonic oscillator potential, is dis- 
cussed in 0. 



* = : E Pn = 



2 



n=l 

00 



I • STn^A - md/2 



71=1 

OO 



i>2:J2 PlK = \{B\k{h bd - e B y- 2 k\B). 



(15) 



n=l 



G. SINGULAR POTENTIALS 

However, if the potential is singular at r — 0, the situa- 
tion changes dramatically. First, only a few first number 
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of the matrix elements (B\Oi\B) are finite. Q 

We can try to proceed along the same lines but then 
have to truncate the series in 1 /5m at the last finite term. 
We want to estimate the difference between this trun- 
cated series and the exclusive sum. 

Let us illustrate this considering a potential with a 
Coulomb behavior at small r, V ~ —a/r, and confining 
at large r. Then (B\k(hbd ~ e B) l k\B) are finite for i < 1, 
but diverge starting from i — 2. We then find that 



r? < 



1 



1 



Such p 2 lead to the estimate 



A 2 / A \ (1+£/a) 



to 2 V 5m 



(17) 



(18) 



More generally, if the above matrix element begins to 
diverge for some value i = K + 1 , the formulas are to be 
replaced by : 



o 2 < 



1 ( 1 



n 1+e \n a 



K+2 



/_A_\ 

m 2 V 5m J 



(K+e/a) 



(19) 



(20) 



Notice that this 5r is smaller than the last retained 



(tt-) 

V dm 1 



K 



term in the OPE series which is of order 
Therefore the 'duality violation' is just smaller than the 
last retained term as for the asymptotic series. This 
means in fact that there is no indication of duality vi- 
olation at this computable order. This is independent of 
a, therefore of the large distance behavior of the potential. 



H. CONCLUSION 

Summarizing our results, the amplitude T(qo, q i ) - the 
T-product, eq. (Q) - can be expanded in inverse pow- 
ers of 5m — q 2 /2m c — qg, the so-called OPE expansion. 
Exact duality would mean that the OPE series was con- 
vergent and equal to T. Actually, this is not exactly the 
case. Even in the favourable case of the regular potentials 
(at r = 0), the OPE series is not convergent, it is only 
asymptotic to the actual T. For singular potentials, the 
coefficients are simply infinite beyond a certain order. 



4 The appearance of infinite coefficients in the OPE series is 
probably due to a breakdown of the power series expansion, 
for instance by fractional powers or logarithms of uiq as seems 
to be the case in the pure Coulomb case |q|. For similar 
phenomena in a perturbation expansion, see |10| . 



Besides these problems concerning the amplitude, ad- 
ditional problems appear for the expansion of the width, 
which is given by a contour integral of T in the q Q com- 
plex plane : the OPE expansion is accurate far from the 
singularities in qa, while the contour has fixed end points 
in the complex plane close to the singularities (Fig. 1). 
In view of this situation, we have computed explicitly the 
difference between the OPE and the actual width. For 
singular potentials, the series must be truncated, and the 
difference is found smaller than the last retained term. 

As to the perspectives opened by this work, we must 
first emphasize that singular potentials seem more inter- 
esting than regular ones. Indeed, in QCD the effective 
quark potential is singular, a smoothed Coulomb singu- 
larity. Moreover, in QCD 2 , one can suspect some sim- 
ilarity with a linear potential \r\, which is also singular 
at the origin in the sense of this paper. For a singular 
potential, we have seen that the entire series must be 
truncated at some order, because the coefficients become 
eventually infinite. We think that such infinite coeffi- 
cients in an entire series expansion correspond to the fact 
that the correct expansion is not entire but must include 
fractional powers and/or logarithms in the expansion pa- 
rameter, i.e. 5m. In QCD, one can argue that the op- 
erator matrix elements are finite due to renormalisation, 
but nevertheless the coefficients still contain logarithms 
of heavy masses. In the non-relativistic case, the object 
of the present paper, the method which has been followed 
does not lead to definite conclusions as regards duality 
for singular potentials : namely, to the order we are able 
to calculate in this paper, we find that there is no dual- 
ity violation, but this leaves open the question of duality 
violation at some higher order Q To proceed further, one 
would have to devise new methods to obtain the above 
conjectured generalized expansions. 
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05 HT 9 HVA3. Laboratoire de Physique Theorique is 
Unite Mixte de Recherche CNRS-UMR8627. 



5 In the context of QCD2, one has demonstrated duality up to 
the order 1 / rriQ and it may be beiieved that duality has been 
fully demonstrated in higher orders j| . However, a comment 
is in order here. In it was shown that the matrix element 
of the leading operator (B\QQ\B) is dual to the sum of the 
widths of the fuli tower of resonances. Therefore, one can 
suspect that there is a difference between the actual width 
and the OPE, that is of higher order 1/mg, corresponding to 
the extrapolated width of the kinematically forbidden states. 
This difference, however, has the same order 1/mo as the 
matrix elements of the higher dimension operators |9| . It was 
then assumed that both quantities are dual to each other, but 
the corresponding OPE coefficients were not calculated and 
we have not found where this assumption was demonstrated. 
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